Finite-temperature phase diagram of the Heisenberg-Kitaev model 
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We discuss the finite-temperature phase diagram of the Heisenberg-Kitaev model on the hexagonal lattice, 
which has been suggested to describe the spin-orbital exchange of the effective spin- 1/2 momenta in the Mott 
insulating Iridate Na2lr03. At zero-temperature this model exhibits magnetically ordered states well beyond 
the isotropic Heisenberg limit as well as an extended gapless spin liquid phase around the highly anisotropic 
Kitaev limit. Using a pseudofermion functional renormalization group (RG) approach we extract both the 
Curie-Weiss scale and the critical ordering scale (for the magnetically ordered states) from the RG flow of the 
magnetic susceptibility The Curie-Weiss scale switches sign - indicating a transition of the dominant exchange 
from antiferromagnetic to ferromagnetic - deep in the magnetically ordered regime. For the latter we find no 
significant frustration, i.e. a substantial suppression of the ordering scale with regard to the Curie-Weiss scale. 
We discuss our results in light of recent experimental susceptibility measurements for Na2lr03. 

PACS numbers: 71.20.Be, 75.25.Dk, 75.30.Et, 75.10.Jm 



In the realm of solid state physics, frustration refers to the 
phenomena that arise from the competition between interac- 
tions that cannot be simultaneously satisfied: typically a large 
degeneracy of ground states and a suppression of thermal or- 
dering by fluctuations 1 1 1. For many magnetic solids a pecu- 
liar form of frustration, so-called geometric frustration, can 
arise when interactions are incompatible with the underlying 
lattice symmetry |2|. A prominent example of the latter are 
spin- 1/2 Heisenberg antiferromagnets on non-bipartite lattice 
structures, for which there is no straight- forward generaliza- 
tion of the Neel state - the common ground state for bipartite 
lattices - but which can instead harbor more exotic ground 
states, including commonly elusive spin liquids |3|. Even 
for bipartite lattices one can encounter geometric frustration 
when considering so-called orbital degrees of freedom, which 
occur in a large class of transition metal oxides that exhibit 
Jahn-Teller ions IH. For the latter crystal field splitting of- 
ten results in a single electron (or hole) occupying the doubly 
degenerate Cg level, for which the orbital occupation is then 
cast in terms of a pseudospin-1/2. In contrast to ordinary spin 
degrees of freedom the exchange interactions between these 
orbital degrees of freedom - arising from Jahn-Teller distor- 
tions and/or superexchange - are highly anisotropic and even 
for simple bipartite lattices cannot be simultaneously satisfied, 
which has been shown to result, e.g. in a non-trivial phase di- 
agram of competing orbital orders on the cubic lattice | 5 1 or 
an orbital Coulomb phase on the diamond lattice 0|. 

In this manuscript, we consider a class of materials, certain 
Iridates, where strong spin-orbit coupling (SOC) results in ef- 
fective degrees of freedom, which fall between the two oppos- 
ing cases above. While Iridates have attracted much recent 
attention as candidate materials for topological insulators Q, 
our study is motivated by a family of materials of the form 
A2lr03, such as Na2lr03, which has recently been shown to 
be a Mott insulator |8|. In these Iridates the Ir^+ (5d^) ions 
form a quasi two-dimensional hexagonal lattice of effective 
j = 1/2 momenta. The latter arise from the combined effect 



of crystal field splitting of the d-orbitals, resulting in a sin- 
gle hole (5 electrons) occupying the lowered t2g orbitals, and 
spin-orbit coupling then giving rise to two Kramers doublets, 
four electrons filling the (lower) j = 3/2 quartet and a single 
electron in the j = 1/2 doublet. The exchange interactions 
between these effective moments have been argued O [TOl to 
reflect both the original spin exchange in terms of an isotropic 
Heisenberg coupling as well as strongly anisotropic orbital in- 
teractions in terms of a Kitaev-type exchange 

Huk[c^] = {I - a)^ai ■ aj - 2a ^ ^i^] ^ (1) 

7 — links 

where the denote the effective spin- 1/2 moment of the Ir^+ 
ions and 7 = x^y^z indicates the three different links of the 
hexagonal lattice. The two couplings are found 1 10] to en- 
ter with opposite sign, i.e. the isotropic exchange is antifer- 
romagnetic, while the anisotropic exchange is ferromagnetic. 
Varying the relative coupling strength < a < 1, the model 
interpolates from the ordinary Heisenberg model with a Neel 
ground state for a = to the Kitaev model for a = 1, which 
even for ferromagnetic interactions is highly frustrated and ex- 
hibits a gapless spin-liquid ground state 1 11 1. One might thus 
wonder how the level of frustration varies between the spin 
and orbital dominated limits of this model. This question is 
also fueled by recent experiments 1 8 1 on Na2lr03 that reported 
magnetic susceptibility measurements, which besides provid- 
ing unambiguous evidence of the effective spin- 1/2 moments 
also reported a considerable suppression for the onset of mag- 
netic correlations below Tn ~ 15 K in comparison with a 
Curie- Weiss temperature of 9cw ^ — 116 K. In particular, 
one might wonder whether this suppression of magnetic or- 
dering might be interpreted as arising from a proximity to the 
highly exotic spin liquid phase of the Kitaev model, despite 
recent resonant x-ray magnetic scattering experiments fT2\ 
reporting indications of a conventionally ordered magnetic 
ground state. 
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FIG. 1. (color online) a) Zero-temperature phase diagram of the Heisenberg-Kitaev model ([T]) exhibiting two magnetically ordered (AFM, 
s-AFM) and a spin liquid phase, b) Positions of the ordering peaks for the magnetically ordered phases in the extended Brillouin zone (BZ). 
The inner hexagon indicates the first BZ. Note that s-AFM order possesses two inequivalent peak positions in the extended BZ. c) Evolution 
of the /c-space resolved static magnetic susceptibility upon variation of the coupling parameter a. 



In this paper, we address the above questions by inves- 
tigating the finite-temperature phase diagram of model ([T]). 
We use a recently developed pseudofermion functional renor- 
malization group (PF-FRG) lfT3llT5ll approach to compute the 
magnetic susceptibility from the pseudofermion two-particle 
vertex function evolving under an RG flow with a frequency 
cutoff A. Neglecting particle fluctuations induced by thermal 
fluctuations, i.e. assuming that we are in the strong coupling 
limit where the on-site repulsion exceeds the temperature, we 
argue and numerically substantiate that the PF-FRG provides 
a suitable tool to obtain both finite-temperature and ground- 
state properties of the model allowing for a direct comparison 
to thermodynamic experiments. In particular, we extract the 
high-temperature Curie-Weiss behavior from the RG flow, the 
onset of magnetic ordering (from the breakdown of the RG 
flow) and momentum-resolved magnetic susceptibility pro- 
files, which also allow to identify the nature of the various 
ground states of model ([T]). 

Numerical simulations - The PF-FRG approach |[T3UT5l 
starts by reformulating the spin Hamiltonian in terms of a 
pseudofermion representation of the spin- 1/2 operators = 
1/2 P =t, I, M = ^, z) with fermionic 

operators and and Pauli-matrices a^. Such a representa- 
tion enables one to apply Wick's theorem leading to standard 
Feynman many -body techniques. In this pseudofermion lan- 
guage, quantum spin models become strong-coupling mod- 
els with zero fermionic bandwidth and a finite interaction 
strength. The major conceptual advancement of the FRG ap- 
proach (T2 IH ri6tflSi is that it allows to tackle this situation 
by providing a systematic scheme for the (infinite) resum- 
mations needed in a direct perturbative attack. These FRG 
summations are obtained in different interaction channels by 
generating equations for the evolution of all m-particle vertex 



functions under the flow of a sharp infrared frequency cutoff 
A. To reduce the infinite hierarchy of equations to a closed 
set, a common approach is to restrict oneself to one-loop di- 
agrams. The PF-FRG extends this approach by also includ- 
ing certain two-loop contributions 1 19] to retain a sufficient 
backfeeding of the self-energy corrections to the two-particle 
vertex evolution. It is this two-particle vertex, which allows 
to compute the magnetic susceptibility - our main diagnostic 
tool to study model ([T]). 

The FRG equations are simultaneously solved on the imag- 
inary frequency axis and in real space. A numerical solution 
requires to i) discretize the frequency dependencies and ii) to 
limit the spatial dependence to a finite cluster, thus keeping 
correlations only up to some maximal length. In our calcu- 
lations the latter typically extends over 7 lattice spacings cor- 
responding to a correlation area (cluster size) of 112 lattice 
sites for the hexagonal lattice at hand. The onset of sponta- 
neous long-range order is signaled by a sudden breakdown of 
the smooth RG flow, while the existence of a stable solution 
indicates the absence of long-range order, see Refs.[T3landfT4l 
for further technical details. 

Zero-temperature states- We start our discussion by first 
considering the zero-temperature phases of the Heisenberg- 
Kitaev model ([T]) and by recapitulating previous T = re- 
sults ifTOl . Interpolating the relative coupling strength a be- 
tween the Heisenberg limit {a = 0) and the Kitaev limit 
(a = 1), a sequence of three phases has been observed 1 10 |: 
The Neel ordered (AFM) state of the Heisenberg limit is stable 
for a < 0.4, when it gives way to a 'stripy' Neel ordered (s- 
AFM) state illustrated in Fig. [T^) which covers the coupling 
regime 0.4 < a < 0.8. In the extended parameter regime 
0.8 < a < 1 the collective ground state is a gapless spin liq- 
uid (SL) where the emerging gapless excitations are Majorana 
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FIG. 2. (color online) The inverse susceptibility 1/x as obtained 
from pseudofermion FRG calculations for various coupling parame- 
ters a. The solid lines indicate fits to a Curie- Weiss law ([2]). 



fermions forming two Dirac cones in momentum space 1 11|. 

We have calculated the characteristic magnetic susceptibil- 
ity profiles for these states within our PF-FRG approach, as 
given for various values of a in the main panel of Fig. [T] 
where we plot the susceptibility just above the (finite) order- 
ing scale Ac below which the RG evolution becomes unsta- 
ble. We adopt the extended Brillouin zone (BZ), illustrated in 
Fig.[T])), appropriate for the two-site unit cell of the hexagonal 
lattice. For the Neel ordered phase we observe characteristic 
corner peaks in the susceptibility. This magnetic AFM sig- 
nature remains robust for the full extent of the phase up to 
a ^ 0.38 where we observe a relatively abrupt shift of the 
susceptibility maxima. The latter is indicative of a first-order 
phase transition, which is in tune with previous T = numer- 
ical studies (TOllJOl. Above a ^ 0.38 we observe the onset of 
the second magnetically ordered phase, the stripy AFM, for 
which the susceptibility signature comes in the form of two 
(not symmetry related) maxima, with a dominant peak along 
/ca^ = in the second BZ and a smaller peak along /c^ = in 
the first BZ. From an extrapolation of the finite-temperature 
crossover line that separates dominant AFM and dominant s- 
AFM fluctuation regimes, we can locate the zero temperature 
transition at a ^ 0.4 in correspondence with previous studies 
CniEOl. Within the s-AFM phase, the point a = 0.5 stands 
out for which the exact quantum ground- state has been shown 
ifTOl to be the classically ordered state (without any dressing). 
In our calculations, the absence of (quantum) fluctuations at 
this point is indicated by remarkably sharp response peaks, 
see Fig.[T]:). For a > 0.5 we observe a pronounced decrease 
of the peaks. In contrast to the transition between the two 
magnetically ordered phases, the phase transition from the s- 
AFM phase to the spin liquid phase around the Kitaev limit is 
more subtle to detect in our calculations. Around the previ- 
ously reported quantum critical point |[T0l[20l at a ^ 0.8 we 
observe a smooth evolution of the susceptibility profile into 
the one expected for the SL phase: a pure cos-type suscepti- 
bility, reminiscent of the purely nearest-neighbor correlations 
in this phase L21J . 



FIG. 3. (color online) The Curie -Weiss scale Acw obtained from 
fitting the inverse susceptibilities in Fig. [2] to a Curie-Weiss law for 
varying coupling parameter a. Around a ^ 0.67 the Curie-Weiss 
scale switches sign indicating a transition of the dominant exchange 
from antiferromagnetic (Acw < 0) to ferromagnetic (Acw > 0). 

Finite -temperature physics - We now turn to the finite- 
temperature properties of model ([T]). To make a connection 
between the flow parameter A and the temperature T, we 
follow a line of thought first discussed by Honerkamp and 
Salmhofer |22|: Both the flow parameter A and the temper- 
ature T did as infrared frequency cutoffs. While the former 
is implemented as a sharp cutoff in the continuous frequency 
space, the latter allows a description in terms of discrete Mat- 
subara frequencies, where the smallest mesh point sets a lower 
bound of the energy resolution. Even though the precise cut- 
off procedures associated with A and T are hence quite dif- 
ferent, we find that the identification of the two scales leads 
to qualitatively correct results; quantitative uncertainties pos- 
sibly enter in our estimates of the ordering instability and its 
critical scale Ac. 

At high temperatures, we find that the homogeneous sus- 
ceptibility calculated from the RG flow for various A nicely 
reproduces the expected Curie-Weiss behavior 

X = C/{A - Acw) , (2) 

as shown in Fig.|2] which allows to extract rather precise nu- 
merical estimates for the Curie- Weiss scale Acw. with the 
latter being plotted in Fig. [3] Notably, we observe that the 
Curie-Weiss scale changes sign - indicating a transition of the 
dominant exchange from antiferromagnetic to ferromagnetic 
- around a ^ 0.68, which is still deep in the magnetically 
ordered regime. Such a change of the dominant exchange is 
already suggestive from a semiclassical analysis of ([T]), which 
gives Tew = — 3/4 + 5a/4 and thus indicates a sign change 
of the Curie- Weiss temperature around a = 0.6. This further 
supports that the cutoff A indeed retains the features of a tem- 
perature parameter and justifies our assumption that A can be 
used to deduce finite-temperature properties of model ([T]). 

We can now return to the question whether a substantial 
frustration builds up as one interpolates between the spin- 
dominated (unfrustrated) Heisenberg regime to the orbital- 
dominated (strongly frustrated) Kitaev regime. A commonly 
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FIG. 4. (color online) Ordering scale Ac obtained from the FRG cal- 
culations for various coupling parameters. The dashed line indicates 
the crossover from dominant AFM to dominant s-AFM fluctuations 
as well as an extrapolation below the ordering transition down to 
T = 0. A regime of enhanced numerical uncertainties is encoun- 
tered near a ^ 0.8. The inset shows the RG flow of the magnetic 
susceptibility versus frequency cutoff A. The arrows indicate the es- 
timated ordering temperatures Ac where the RG flow breaks down. 

used measure for frustration is the ratio between Curie- Weiss 
and ordering scale, the so-called frustration parameter 

/= lOcwI/Te^ |Acw|/Ac, (3) 

with a small value / ^ 5 indicating the absence of frustra- 
tion, and systems with / ^ 10 being commonly referred to as 
highly frustrated 1 1 1. We estimate the ordering scale Ac from 
the breakdown of the RG flow, as shown in the inset of Fig.|4j 
with the result being plotted in the main panel of Fig.[4]for the 
full parameter range a except for a region around the transi- 
tion between the s-AFM phase and the SL (a ^ 0.8), where 
our approach does not allow to reliably calculate the transition 
temperature. As shown in Fig. [5] we observe a rather constant 
plateau / ^ 2 for the frustration parameter in the AFM regime 
before it decreases linearly starting around a 0.4 with the 
onset of the s-AFM phase. At a = 0.5 in the s-AFM phase, 
model ([T]) can be mapped to a fluctuation-free classical system 
ifTOl . This is consistent with our result of/^lata = 0.5, 
as a frustration parameter close to unity signals the absence 
of fluctuation-induced frustration. For larger a the frustration 
parameter goes through zero as the Curie- Weiss scale changes 
sign and beyond a regime of numerical uncertainty rapidly di- 
verges as expected for the spin-liquid regime, see Fig. [5] 

Connection to experiments - We close our discussion of 
the finite-temperature properties of model ([T]) by comparing 
our findings to recent thermodynamic measurements [S] on 
the Iridate Na2lr03. The reported Curie- Weiss temperature 
of 6cw ^ — 116 K indicates a dominant AFM exchange and 
the considerable suppression of magnetic ordering down to 
Tn ~ 15 K corresponds to a frustration of / ^ 8. While our 
finite-temperature analysis of the Heisenberg-Kitaev model 
([T]) indicates an AFM Curie- Weiss temperature for a wide 
range of couplings < a < 0.68, the ground states in this 
regime are relatively simple, magnetically ordered states that 



FIG. 5. (color online) The frustration parameter, i.e. the ratio of 
Curie-Weiss temperature Acw and ordering temperature Ac. 

do not give rise to a significant suppression of the ordering 
temperature with the frustration parameter / never exceeding 
/ ^ 2 in our calculations for this regime, see Fig. |5] On the 
other hand, we find a strong suppression of the ordering tem- 
perature in the spin-liquid phase for a > 0.8 and its proximity, 
but the dominant couplings in this regime are ferromagnetic 
(6cw > 0). To reconcile the combination of an AFM Curie- 
Weiss temperature and a simultaneous suppression of the or- 
dering temperature, one might thus want to look beyond the 
Heisenberg-Kitaev model ([T]). In particular, one might want 
to consider various mechanisms that could suppress the order- 
ing temperature in the magnetically ordered regime, such as 
a next-nearest neighbor exchange introducing geometric frus- 
tration as suggested in Ref. 23 or the role of disorder |24|, es- 
pecially in the form of non-magnetic impurities arising from 
the experimentally observed 1 8 1 site mixing between Ir and Na 
atoms. While the current analysis might suggest that Na2lr03 
is not in close proximity to the spin liquid phase of the Kitaev 
limit and its cousin - a topological spin liquid in the pres- 
ence of a magnetic field |20|, one might still speculate how 
one could drive the system closer to that regime. One promis- 
ing path to experimentally increase the anisotropic couplings 
might be to exert pressure along the a6-plane to counteract the 
c-axis lattice distortion in the material, which quenches the 
SOC. A similar relief of the lattice distortions might also be 
expected when replacing the Na ions by smaller Li ions ll25l 
and consider Li2lr03 as a candidate material for more exotic 
ground states. 
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